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ABSTRACT 

The MBH~Afhost relation in quasars has been probed only in a limited parameter space, namely at 
Mbk^ 10^ Mq and Mhost~ 10^^ M©. Here we present a study of 26 quasars laying in the low-mass 
end of the relation, down to Mbh^ lO"^ Mq. We selected quasars from the SDSS and HST-FOS 
archives, requiring modest Mbh (as derived through the virial paradigm). We imaged our sources in 
H band from the Nordic Optical Telescope. The quasar host galaxies have been resolved in 25 out 
of 26 observed targets. Host galaxy luminosities and stellar masses are computed, under reasonable 
assumptions on their star formation histories. Combining these results with those from our previous 
studies, we manage to extend the sampled parameter space of the MeH-Afhost relation in quasars. 
The relation holds over 2 dex in both the parameters, similarly to what observed in low-luminosity 
AGN and in quiescent galaxies. For the first time, we are able to measure the slope of the MeH-Afhost 
relation in quasars. We find that it is consistent with the linear case (similarly to what observed in 
quiescent galaxies). We do not find any evidence of a population of massive black holes lying below 
the relation. 

Subject headings: quasars: 



1. INTRODUCTION 



Massive black holes (BHs) ar e ubiquitously found in 
the c e ntre of massive galaxies (jKormendv fc Richstond 
[Tool IDecarh et all [20071: IGallo et all \^oW . Their 
masses (Mbh) show strong correlations with large-scale 
properties of the host galaxi es, namely the stel 



lar velocity dispe rsi on (a 
Ferrarese fc Merrittl 120001 : 



IGebhardt et al.l 
iGultekin et al.l 



2000; 



200S; 



Graham et al.l 12011 ). the luminosity and mass of the 



spheroidal compone nt (£host i -Mhost; s ee'Mago rrian et al.l 
[1998; Mar coni fc Hu nt 2003: IHaringfc Rix 20Q1). These 
relations have been interpreted as the outcome of a 
joint evolution between B Hs and their host g alaxies 
JSilk fc Rccs 1998; Kauff mann fc Haeneltl 12000: 
2Q05i: iWvithe fc Loeb ,20061: iRober 
Hopkins et al.l l2007t iMalbon et all 



Kine 

Robertson et al.l |200€ 



2007r 




20091) : In this scenario, the growth of BHs through 



accretion regulates the gas cooling in the outskirt of 
host galaxies through energy or momentum injection 
(feedback), thus quenching the formation of stars. 
Galaxy mergers may also play a role in this scenario, 
as gravitationally-induced dynamical instabilities may 
trigger both star format ion bursts and gas inflows 
fuelling the BH activity ([K auffmann fc Hacnclt 2000'; 
Di Matteo. Soringel fc Hernauist .2005.; Canalizo et al.. 



2007t iBennert et al.ll2008l although see also Cisternas et 



D- 



al. 2011) 
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The MeH-host galaxy relations have been pinned down 
on an albeit small set of local, mostly quiescent galax- 
ies. The sampled parameter space ranges over 3 dex in 
masses, from a few million to a few billion solar masses 
in terms of Mbh- Extending these studies beyond the 
local Universe is challenging. On one side, the influence 
radius of BHs, i?inf, i.e., the radius where the gravita- 
tional potential is dominated by the singularity, is re- 
solved only in very nearby objects (distances < few tens 
Mpc) with high Mbh values. For any other sources, indi- 
rect tracers of Mbh are required. The most commonly- 
adopted indirect estimator of Mbh is based on the width 
of broad emission lines and the size of the broad line 
region (BLR). This can be done o nly in type-1 AGN, 
where broad lines are observed (IPeterson fc Wandell 
2000; Vcstcrgaard |00l iVestergaard fc PetersonI " [2006r 
Decarli et al.. 2010d^ . This approach allows to esti- 
mate Mbh from single-epoch spectra in ^ 100 OOP 
quasars up to 2 ~ 5 from SDSS spectra (|Shen et al.1 
2010), and in most of the z ~ 6 quasars kno wn to 
day (|Willott McLure fc .Tarvisl [2 003*: K urk et al.l [20071 : 
iWillott et aL[[2OT0l: iDe Rosa et al.. 2011. H 

On the other side, the properties of host galaxies 
are hard to measure in distant sources. Bright active 
nuclei (necessary to measure Mbh) can easily outshine 
the light of their host galaxies . Observations i n ex- 
cellent seeing cond itions (e.g., iKotilainen et al.l 120091 : 
Targett et al.l 12011 ) or based on adapti ve optics (e.g. 



Groom et alTl2004t iFalomo et~aL[ 120051 ) are required. 



Diffraction-limited observations with HST have also 



^ Caveats to this technique aris e as BLR clouds may 
be supported by radiation pressure lIMarconi et al.l 120091 ) or 
move non-virially; projection effects depending on the geom- 
etry and orientation of the BLR may hinder our abi lity to 
actual l y measure the orbital velocit y of clouds I Deca rli et al.l 
I2008bl : IDecarli, Dotti fc Treve^ [20TlD : different emission lines 
may be produced in regions where t h e gas d ynamics are differ- 
ent IINetzer e t al.' '200?; Decarli ct al. 2008a; 'L abita et aP l2009l : 
IFine. Jarvis fc Mauch 2011: Richards ct al. 2011 ). 
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signif ica ntly contribu te d in th is field (iBahcall et al] 



1997t iKukula et al I [200lt iRidewav et all 12001 



Hamilton et all 120021: iDunlop et al l i2003t iFlovd et al 



2004; Pe ng et alJl2q06at iKim et al.M2q08HJalinke et al 



2009t IBennert et alJ 120111 : I Jiang et all 1201 It ), although 



some concerns a bout the reprod ucibility of the PSF 
have been arised (|Kim "eFalllSool . Up to now, ~ 300 
quasar host galaxies have been resolved up to z » 3 , 
and most of them at z < 0.5 (see iKotilainen et al.ll2009l 
and references therein). 

In order to understand the processes and timescales 
which led to the onset of the BH-host galaxy relations, 
two key observational tests are required. The first one 
consists in tracing the evolution of the BH-to-host galaxy 
mass ratio (T = MBH/^^hostJ as a function of Cosmic 
Time. If, e.g., T{z > 0) < T{z « 0), then we can argue 
that the BHs in these systems still have to accrete in al- 
ready formed bulges. Vice versa, T{z > 0) > r(z « 0) 
could suggest a rapid growth of the BHs, followed by 
a slower build-up of the spheroids. Most of the studies 
on the evolution of the BH-host galaxy relations sug- 
gest that at high redshift, for a given mass of the host 
spheroid, th e harbore d BH is more massiv e than at low-z 
(iMcLiire et al. 2006t iPeng e t al. 2006a b !: iMerloni et al.l 
iDecarli et alTl2010bl: ' iBennert et all 1201 1[ ). with 
cx (1 -I- z)'^/'^. It is interesting to note that the 
host galaxy of J1148-f5251, the highest-^ SDSS quasar, 
z = 6.42, sho ws a A/BH~host galaxy dy namical mass 



[2oTa 

T{z) 



at 

ratio of 



zi, sno ws a iviBH~riost galaxy gy i 
0.13 (1 Walter et al.ll2003l l200l . 



in an order- 



of-magnitude agreement with the extrapolation from the 
z < 3 studies. 

A second observational test to probe the onset of the 
BH-host galaxy relations is to trace the low-mass end 
of the BH-host galaxy relations. Different initial host 
galaxy mass, BH seed mass and build-up processes pro- 
duce different slopes of the relations, especially at low- 
mass. Light seeds (10 — 100 M©) are expected as the 
remnants of metal-pure stars in the early Universe, while 
heavier seeds (up to ^ 10^ M©) can result from the direct 
collapse of primordial gas clouds. The former ones would 
produce a larger scatter in the MBH~liost relations, a 
higher occupation fration in relatively small galaxies and 
a lower cutoff in the minimum Mbh with res pect to the 
latter (see, e.g.. iVolonteri fc Nataraianll2009l and refer- 
ences therein). Some authors even claim that the Mbh~ 
(7, relation itself is just the upper limit of a broader 
distribution, with a number of (hard-to-detect) modest- 
mass BHs embed ded in relatively massive galaxies (e.g. 
IBatcheldod[20l0l) . 

The low-mass end of the MBH-host galaxy relations 
has been probed down to few 10^ Mbh in quiescent galax- 
iesor in low-luminosity AGN at low-z fiGreene fc Hoi 
200l [2007t IDongetall [20071: iTh orton et all I2OOI 



Xiao et al. [120111 : I Jiang et al. suggesting that the 



MBH-host relations hold in quiescent or mildly active 
galaxies. However, no effort has been attempted so far to 
extend this test to higher luminosity AGN. Quasars are 
ideal probes of the BH~host galaxy relations at z > 0. 
However, while F in quasars has been measured up to 
very high redshift, the ranges of Mbh and Mhost inves- 
tigated up to date are limited, and comparable with the 
observed scatter of the BH-host galaxy relations. Filling 
the low-mass end of the MBH-Mhost relation therefore 
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Fig. 1. — Comparison among the Mbh distributions in the SDSS 
quasars at 2 < 0.5 (Shcn et al. 2010, ; top panel), in our 2010 study 
{middle panel) and in the present work {bottom panel). The latter 
is clearly more representative of the general quasar population at 
low redshift. 



represents a main challenge and a fundamental step in 
our comprehension of the BH-host galaxy evolution. In 
this paper, we present ground-based NIR observations of 
quasars at z < 0.5 selected so that virial Mbh < 10^ 
A/bh- Our imaging campaign successfully resolved 25 
quasar host galaxies. This enables us to directly probe 
the slope of the Msn-Mhost relation in quasars. 

The structure of this work is the following: in Section[2] 
we describe the sample. In Section|3]we present the anal- 
ysis of the spectra and we derive Mbh in all our sources. 
The new observations, the data reduction and analysis 
and the results from the imaging campaign are presented 
in Section [D In Section [5] we discuss our results. Con- 
clusions are summarized in Section [51 Throughout the 
paper we will assume a standard cosmology with Hq — 70 
km s~i Mpc^i, f^in = 0.3 and Qa = 0.7. 

2. THE SAMPLE 

We selected quasars at z < 0.5 with availabl e H/3 obser- 
vatio ns in the Sloan Digital Sky Survey (SDSS l^rk et al.l 
l2000f ) spectroscopic database or CiV or Mgii obser- 
vations in the HST-Faint Object Spectrograph (FOS) 
archive. We require that black hole virial masses, com- 
puted as described in Section [3j range between 10 and 
10^ Mq. Note that, out of t he 62 z < 0.5 quasa rs ex- 
amined in our previous study (|Decarli et al.ll2010al ) , only 
8 (13%) had Mbh< 10^ Mq. We then selected quasars 
having at least 3 relatively bright (m/f =10-15 mag) stars 
within 2' (corresponding to the half size of the NOT- 
Cam field of view), and a number of fainter stars in or- 
der to have an accurate description of the Point Spread 
Function (PSF). Observability constraints and modest 
weather losses further limited our analysis to 26 targets 
(see Table [1]) , mostly detected in radio wavelengths ac- 
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TABLE 1 

The sample. (1) Target name. (2-3) Right ascension and 
dechnation (J2000). (4) Is the target detected in radio 
wavelengths, according to the Vcron-Cctty & Vcron ( 2010) 
catalogue? (5) Catalogue redshift. (6) Apparent visual magnitude 
of the quasar. 



Name 



(1) 



RA 



(2) 



DEC 



(3) 



Radio? 



(4) 



(5) 



V 

[mag] 
(6) 



PB5723 


00 


05 


47.5 


+02 


03 


02 


N 


0.234 


16.60 


PG0026+12 


00 


29 


13.7 


+ 13 


16 


04 


Y 


0.145 


15.41 


PG0052+251 


00 


54 


52.1 


+25 


25 


39 


Y 


0.155 


15.43 


B20110+29 


01 


13 


24.2 


+29 


58 


16 


Y 


0.363 


17.00 


PKS0214+10 


02 


17 


07.6 


+11 


04 


10 


Y 


0.408 


16.46 


J02321+0008 


02 


32 


11.8 


+00 


08 


03 


Y 


0.432 


19.10 


J02331-0909 


02 


33 


10.6 


-09 


09 


40 


Y 


0.388 


18.45 


J03010+0004 


03 


01 


00.2 


+00 


04 


29 


Y 


0.486 


19.33 


J03323+0106 


03 


32 


18.0 


+01 


06 


48 


N 


0.482 


18.91 


J03579-0550 


03 


57 


59.0 


-05 


50 


15 


Y 


0.439 


18.93 


B20752-I-25A 


07 


55 


37.0 


+25 


42 


39 


Y 


0.446 


18.00 


J08044+1904 


08 


04 


42.1 


+19 


04 


26 


Y 


0.346 


19.40 


J08285+2748 


08 


28 


53.5 


+27 


48 


33 


N 


0.330 


20.00 


J08305+0802 


08 


30 


57.4 


+08 


02 


34 


Y 


0.319 


19.20 


PG0844-I-349 


08 


47 


42.4 


+34 


45 


03 


Y 


0.064 


14.50 


J09010+3538 


09 


01 


00.9 


+35 


38 


09 


N 


0.302 


19.20 


PG0947+396 


09 


50 


48.3 


+39 


26 


51 


Y 


0.206 


16.39 


PG0953+415 


09 


56 


52.4 


+41 


15 


23 


Y 


0.234 


15.32 


TON1187 


10 


13 


03.1 


+35 


51 


23 


N 


0.079 


14.75 


TEX1156+213 


11 


59 


26.2 


+21 


06 


56 


Y 


0.349 


16.90 


Q1214+1804 


12 


16 


49.0 


+ 17 


48 


04 


N 


0.374 


17.30 


PG1404-I-226 


14 


06 


21.9 


+22 


23 


47 


Y 


0.098 


15.82 


PG1415+451 


14 


17 


00.8 


+44 


56 


06 


Y 


0.114 


15.24 


PG1626+554 


16 


27 


56.1 


+55 


22 


31 


Y 


0.132 


15.68 


4C73.18 


19 


27 


48.5 


+73 


58 


02 


Y 


0.302 


16.50 


PKS2251+11 


22 


54 


10.4 


+11 


36 


39 


Y 


0.325 


15.82 



cording to the iVeron-Cettv fc VeronI ()2010D catalogue. 
This new sample is then matched with the 62 z < 0.5 
targets present ed in our pr evious study of the Mbh- 
Mhost relation (jPecarli et a l. 2010a, b). Figure [T] com- 
pares the distribution of A/bh in the general SDSS sam- 
ple at z < 0.5, computed from the conti nuum lumi- 
nositie s and H/3 width estimates r eported by iShen et alj 
(|2010t) . with the ones presented in iDecarli et al.l ()2010aD 
and in the present study. 

3. THE SPECTROSCOPIC DATASET 

The spectroscopic dataset consists of pipeline- 
processed, publicly available spectra from the SDSS or 
FOS archive. SDSS spectra have A/AA 2000 and a 
spectral range between 3800 and 9000 A. Uncertainties 
on wavelength calibration amount to 0.05 A, while flux 
calibration formal errors account to 5%. FOS spectra are 
taken from the compilation of re-calibrated quasar and 
AGN FOS spectra by Evans & Koratkar (2004). Obser- 
vations were performed with a number of different grat- 
ings both at high spectral resolution (1-6 A diode" ^, 
A/AA ~ 1300) and low spectral resolution (6-25 A 
diode"^, A/AA w 250) covering various spectral win- 
dows from 1140 A to 3275 A. Photometric uncertain- 
ties are usually 5-10%, while typical wavelength cali- 
bration un certainties are around 0.5 channels rms (see 
lEvans &: Koratkar.,2004 ). 
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Fig. 2. — Examples of H/3 line fitting. The observed spectra, 
shifted to the rest frame, are plotted in solid, black lines. The 
various components of the models are shown: the power-law (black, 
dotted line); the host galaxy and Fell templates (green, solid line); 
the line model (red, solid line) and its components (red, dotted 
lines). Fit residuals are shown in blue, dashed lines. 



Spectral analysis follo w s the same approac h de- 
scribe d in 'Deca rli et al.l ()2010a[ ) and ^Dc Ros a et al.l 
(j2011| ). Briefly, continua are modeled as a superposi- 
tion of a power-law component, the host galaxy emis- 
sio n (only in the optical spectra; the Elliptical model 
by iMannucci et al.l I2001L was adopted as a template) 
and the Fell multiplets (m odeled on the template by 
IVestergaard fc WilkesI [200l . Relevant broad emission 
lines (H/3, Mg li and C iv) are modeled with the super- 
position of two gaussian curves with the same peak wave- 
length. The spectral resolution of the SDSS data is ad- 
equate to allow an easy identification of the narrow H/3 
emission with respect to the broad line. Examples of the 
line fitting are shown in Figure [2] 

We measure the continuum luminosity at 1350, 3000 
and 5100 A from the fitted power-law. The broad line 
luminosity and FWHM a re measured from the line model 
(see iDecarli et al.l l2008a[ for a discussion on the fitting 
technique and the parametrization of the relevant quan- 
tities). Virial black hole masses are comp uted with the 
same recipies used in IDecarli et al.l (|2010al) . adopting ge- 
ometrical factors of 1.6 for H/3 and Mg ii and 2.4 for C IV. 
Table [2] lists the main measurements and inferred quan- 
tities from the spectroscopic analysisQ- 

4. THE IMAGING DATASET 

All the objects in our study have been observed in H 
band in a campaign at the 2.5m Nordic Optical Telescope 

^ We note that, after a careful re-analysis of the spectra, two ob- 
jects (B20110+29 and B20752+25A) show values of A/bh slightly 



exceeding the initial selection criteria, 
these sources in the present analysis. 



Nevertheless, we include 
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Fig. 3. — An example of the analysis of the quasar host galaxies, shown in the case of the quasar J08285+2748. Top panels: Observed 
quasar image (left), contour plot (center) and light profile (right; the PSF and galaxy models are plotted in dotted and dashed lines 
respectively). Central panels: similarly, the images (left), contour plots (center) and light profile (right) of the PSF and the galaxy model. 
Bottom panels: image (left), contour plot (center) and light profile (top right) of residuals after model subtraction. The radial distribution 
of values is also plotted (bottom right). 
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TABLE 2 

Results from the spectroscopic analysis. (1) Quasar name. (2) 
Redshift. (3) Line used in the A/bh estimate. (4) Continuum 
monochromatic luminosity at 1350 (for Civ), 3000 (for Mgll) or 
5100 A (for H^), in erg s"!. (5) Line FWHM in km s^^. (6) 
Virial estimate of the black hole mass, in solar units. 



Name 


z 


Line 


log \Lx 


FWHM 


log A/bh 








[erg s"-*-] 


[km s-i] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


PB5723 


0.234 


Civ 


44.71 


3715 


8.15 


PG0026+12 


0.145 


Civ 


45.22 


2062 


7.92 


PG0052+251 


0.155 


Civ 


45.33 


5914 


8.90 


B20110+29 


0.363 


H/3 


44.81 


6149 


9.33 


PKS0214+10 


0.408 


Civ 


45.71 


4122 


8.79 


J02321+0008 


0.432 


H/3 


44.41 


1727 


7.56 


J02331-0909 


0.388 


H/3 


44.63 


1863 


7.77 


J03010+0004 


0.486 


H/3 


44.45 


6634 


8.76 


J03323+0106 


0.482 


H/3 


44.77 


4282 


8.59 


J03579-0550 


0.439 


H/3 


44.62 


4005 


8.43 


B20752+25A 


0.446 


H/3 


45.06 


9738 


9.50 


J08044+1904 


0.346 


H/3 


44.09 


7322 


8.60 


J08285+2748 


0.330 


H/3 


44.14 


5385 


8.37 


J08305+0802 


0.319 


H/3 


44.14 


6149 


8.48 


PG0844+349 


0.064 


Mgii 


44.57 


3209 


8.20 


J09010+3538 


0.302 


H/3 


44.29 


8495 


8.86 


PG0947+396 


0.206 


Civ 


45.17 


4090 


8.49 


PG0953+415 


0.234 


Civ 


45.45 


3490 


8.50 


TON1187 


0.079 


H/3 


44.20 


2141 


7.60 


TEX1156+213 


0.349 


H/3 


44.93 


5663 


8.94 


Q1214+1804 


0.374 


H/3 


45.04 


3728 


8.65 


PG1404+226 


0.098 


H/3 


44.14 


1036 


6.93 


PG1415+451 


0.114 


H/3 


44.19 


3244 


7.96 


PG1626+554 


0.132 


Civ 


45.18 


4057 


8.49 


4C73.18 


0.302 


Civ 


45.92 


4155 


8.91 


PKS2251+11 


0.325 


Civ 


45.54 


5028 


8.87 



(Roque de Los Muchachos, Spain) using NOTCam. Ob- 
servations were carried out during three observing runs 
in May 2007, April and November 2008. The average 
seeing in B band was 0.7". The 1024 x 1024 pixel NOT- 
Cam detector has a pixel scale of 0.235" pxl~^, yielding 
a field of view size of ~ 4 x 4 arcmin^. Usual observ- 
ing techniques for broad-band NIR imaging of point-like 
sources were adopted. Observations were split in 1 min 
long individual frames. Random jittering patterns within 
a 20" box were adopted in order to perform optimal sky 
subtraction. The total time on each source was 45 min. 

Data reduction was performed using IRAP0. Bad pix- 
els were corrected for in each image using a mask made 
from the ratio of two sky fiats with different illumina- 
tion levels. Sky subtraction was performed for each sci- 
ence image using a median averaged frame of all the 
other temporally close frames in a grid of eight expo- 
sures. Flat fielding was made using normalized median 
averaged twilight sky frames with different illumination 
levels. Finally, images were aligned to sub-pixel accu- 
racy using field stars as reference points and combined 
after removing spurious pixel values to obtain the final re- 
duced co-added image. Zero point calibration is achieved 

* IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



by cross-matching the photometry of field stars with the 
2MASS database. This approach yields typical photo- 
metric uncertainties <0.1 mag. 

We analyzed our data using the Astronomical 
Image Decompositi on an d Analysis package (AIDA; 
[Uslcn ghi fc Falomg |2008[) . an IDL-based software de- 
signed to perform two-dimensional model fitting of 
q uasar images. De t ails on the procedure are presented 
in lKotilainen et ahl ()2007l [20091 ) . and briefly summarized 
here. 

A careful modeling of the Point Spread Function (PSF) 
is crucial to disentangle the extended host galaxy light 
from the nuclear source. To model the PSF shape, we 
used suitable field stars. Each star was modeled with 
four two-dimensional Gaussians, representing the core of 
the PSF, and an exponential feature, representing the ex- 
tended wing of the PSF. Regions contaminated by nearby 
sources, saturated pixels and other defects affecting the 
images were masked out. 

In order to discriminate between resolved and unre- 
solved targets, we first fit the images of our sources with 
the pure PSF model. In most of the cases, an extended 
halo was clearly observed in the residuals. We then re- 
performed the fits using a Scrsic law, (describing the host 
galaxy) plus a point-source (the nucleus), convolved to 
the PSF model. In all but one case (J02331-0909), this 
second fit is significantly better than the fit with the pure 
PSF, as confirmed by the ratio between the two fits 
(see table [3|). 

An example of the outputs of our analysis is shown in 
Figured 

Host galaxy apparent -ff-band magnitudes are then 
converted into rest- f rame .R-band. We use the Elliptical 
galaxy template by iMannucci et al.l (|2001| ) to estimate 
fc-corrections. The host galaxy i?-band absolute magni- 
tude is then converted into a stellar mass by adopting 
the mass-to-light ratio [M / L) of a single stellar popula- 
tion originate d at ^burst = 5 and pa ssively evolving down 
to z = fsee iDecarli et al.l l2010b'. for details). Table [3] 
summarizes the relevant results from the modeling of the 
quasar host galaxies described here. 



5. DISCUSSION 

In FigurelHwe show the MBH--^host relation for quasar 
host galaxies at z < 0.5. The dataset (62 objects from 
Decarli et al. 2010, plus 25 objects with resolved host 
galaxies from the present study) span over 2 dex both 
in A/bh and Mhost- The same MBH~-^host relation ob- 
served for inactive galaxies appears to hold through all 
the sampled range, from 3 x 10^0 to 3 x 10 Mq 
in terms of Mhost- We find that (logF) = -2.843 (in 
excellent ag reement with the Mbh =0-0015 Mhost value 
reported bv lMarconi fc Huntll2003l . for inactive galaxies) 
with a 0.44 dex scatter. Only 3 sources (J023214-0008, 
PG1404-j-226 from this study; lOOH-291 from the old 
sample) lie more than 2-a below the relation. Since the 
sampled parameter space is about 5 times larger than 
the dispersion of the relation, we can exclude that the 
observed MBH-Mhost relation is the upper envelope of a 
population of quasars with relatively small black holes 
hosted by very massive galaxies. The best bilinear re- 
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TABLE 3 

Results from the imaging analysis. (1) Quasar name. (2) Redshift. (3) Apparent observed H-band magnitude of the nucleus. (4) 
Apparent observed H-hand magnitude of the host galaxy. (5) Apparent central surface brightness, as extrapolated from the host galaxy 
model. (6) Chi square ratio between the best fit using a pure PSF and the best fit using a galaxy+quasar model. (7) Effective radius of 
the host galaxy. (8) EUipticity of the host galaxy model. (9) Sersic index of the host galaxy model. (10) Rest-frame R - observed H color 
correction. (11) Resulting absolute rest-frame /?-band magnitude of the host galaxy. (12) Adopted mass-to-light ratio. (13) Stellar mass 

of the host galaxy. (14) BH-to-host galaxy mass ratio. 



Name 


z 


^nuc 




MO 


Xpsf/Xgq 


Re 


Ell. 




R-H 


Mr 


log M/L 


log Mhost 


logP 






[mag] 


[mag] 


[mag "-2] 






["] 




[mag] 


[mag] 


[Mq/Lq] 


[Mq] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11 


(12) 


(13) 


(14) 


PB5723 


0.234 


15.99 


15.81 


5.41 


1.41 


0.32 


0.57 


5.00 


2.59 


-21.93 


0.71 


11.38 


-3.61 


PG0026+12 


0.145 


12.90 


14.94 


16.10 


1.77 


1.44 


0.15 


1.20 


2.58 


-21.64 


0.75 


11.30 


-3.76 


PG0052-f251 


0.155 


13.46 


14.37 


16.34 


6.55 


2.19 


0.19 


1.22 


2.58 


-22.38 


0.74 


11.59 


-3.07 


B20110-I-29 


0.363 


16.04 


16.21 


12.07 


2.51 


3.89 


0.09 


5.00 


2.60 


-22.61 


0.67 


11.61 


-2.48 


PKS0214-flO 


0.408 


14.60 


16.10 


13.98 


1.88 


1.16 


0.00 


2.77 


2.59 


-23.03 


0.65 


11.77 


-3.36 


J02321+0008 


0.432 


99.90 


16.80 


9.93 


1.43 


0.37 


0.00 


3.89 


2.58 


-22.50 


0.64 


11.54 


-4.18 


J02331-0909 


0.388 


15.98 


- 


- 


1.01 


- 


- 


- 


- 


- 


- 


- 


- 


J03010-I-0004 


0.486 


99.90 


17.35 


13.42 


1.17 


0.26 


0.24 


1.89 


2.57 


-22.25 


0.63 


11.43 


-2.87 


J03323-I-0106 


0.482 


18.66 


17.15 


11.55 


1.91 


0.54 


0.61 


3.17 


2.57 


-22.43 


0.63 


11.50 


-3.11 


J03579-0550 


0.439 


18.78 


17.22 


8.06 


1.14 


0.45 


0.35 


5.00 


2.58 


-22.12 


0.64 


11.39 


-3.16 


B20752+25A 


0.446 


14.80 


16.14 


12.51 


1.39 


0.71 


0.07 


2.95 


2.58 


-23.24 


0.64 


11.83 


-2.53 


J08044-I-1904 


0.346 


17.07 


17.40 


17.73 


1.55 


0.98 


0.37 


1.00 


2.60 


-21.30 


0.67 


11.09 


-2.69 


J08285-I-2748 


0.330 


16.98 


17.20 


16.63 


1.58 


0.81 


0.06 


1.53 


2.60 


-21.38 


0.68 


11.13 


-2.96 


J08305-I-0802 


0.319 


17.50 


18.25 


17.66 


1.19 


0.57 


0.33 


0.90 


2.60 


-20.25 


0.68 


10.68 


-2.40 


PG0844+349 


0.064 


13.34 


14.59 


17.55 


1.92 


2.82 


0.27 


0.90 


2.56 


-20.13 


0.78 


10.73 


-2.73 


J09010-I-3538 


0.302 


17.08 


16.35 


14.10 


2.79 


0.91 


0.20 


2.43 


2.60 


-22.01 


0.69 


11.39 


-2.73 


PG0947+396 


0.206 


14.42 


15.18 


16.17 


4.45 


1.92 


0.20 


1.59 


2.59 


-22.24 


0.72 


11.52 


-3.41 


PG0953+415 


0.234 


12.93 


15.80 


19.36 


1.33 


3.39 


0.12 


0.90 


2.59 


-21.93 


0.71 


11.38 


-3.26 


TON1187 


0.079 


13.97 


14.95 


16.83 


3.70 


1.91 


0.20 


1.10 


2.56 


-20.25 


0.77 


10.77 


-3.37 


TEXl 156-1-213 


0.349 


15.38 


15.78 


14.86 


2.50 


1.07 


0.12 


1.98 


2.60 


-22.95 


0.67 


11.75 


-3.01 


Q1214-I-1804 


0.374 


16.36 


16.48 


12.05 


1.11 


0.25 


0.40 


1.97 


2.59 


-22.44 


0.66 


11.54 


-3.09 


PG1404+226 


0.098 


14.35 


14.70 


17.08 


7.62 


2.99 


0.52 


1.04 


2.57 


-20.97 


0.77 


11.05 


-4.32 


PG1415+451 


0.114 


13.65 


14.02 


15.22 


10.56 


1.50 


0.04 


1.30 


2.57 


-22.01 


0.76 


11.46 


-3.70 


PG1626-f554 


0.132 


13.80 


14.78 


16.23 


8.55 


2.04 


0.04 


1.53 


2.57 


-21.60 


0.75 


11.29 


-3.18 


4C73.18 


0.302 


13.55 


16.42 


17.58 


1.44 


1.16 


0.18 


0.90 


2.60 


-21.94 


0.69 


11.36 


-2.83 


PKS2251+11 


0.325 


13.53 


15.70 


18.62 


2.94 


2.43 


0.05 


0.90 


2.60 


-22.86 


0.68 


11.72 


-3.23 



gression fit of the relation is: 



as follows: 



log 



M, 



BH 



M, 



109 ]VI, 







(1.26±0.29)xlog ^ ^ -f (0.04±0.03) 

(1) 

consistent with the relations with a constant Mbh / A/host 
ratio, as observed in quiescent galaxies in the loc al Uni- 
verse ijMarconi fc Huntl[2?M IHaring fc Rb3 [200l. 

When considering subsets of our data, the high-mass 
end shows a slightly smaller scatter (see Table |4]). 

For twelve objects the host galaxies are found to be 
best described by a Sersic law with small index (rig < 
1.5), suggesting the presence of significant disc compo- 
nents. In particular, PG1404-I-226 (incidentally, the ob- 
ject showing the smallest Mbh and the largest deviation 
with respect to the MeH-Mhost relation in our sample) 
shows clear spiral arms in the residuals of the host galaxy 
model. From local galaxy studies, Mbh is found to be 
more sensitive to the properties of the spheroidal stellar 
component rather than of the whole galaxy. On the other 
hand, a bulge+disc decomposition is practically impos- 
sible with ground-based images of quasar host galaxies 
at relatively high redshift. Here we attempt a rule-of- 
thumb correction starting from the Sersic index value. 
We assume that the bulge-to-total luminosity ratio in 
the rest-frame R band, B/T, scales with the Sersic index 



B/T = 



(us - 0.5)/3.5 if < 4 
1 if > 4 



(2) 



This simple analytical form roughly traces the bulk of th e 
B/T values for Us < 4, as found bv lSimard et all ()201lD . 
who performed accurate image deconvolution for ~ 1 mil- 
lion galaxies from the SDSS. Furthermore, it is consistent 
with the operative hypothesis that all the galaxies well 
described by a de Vaucouleurs profile = 4) are bulge 
dominated {B/T w 1), as assumed in our previous study. 
The effect of this correction is to move disc-dominated 
host galaxies towards the left side of Figure ID In par- 
ticular, all but one source at Mbuigo < 10^° Mq would 
lie above the local relation. The best fit relation is in- 
deed flatter (0.88 ± 0.18 instead of 1.26 ± 0.29), but still 
consistent with the linear case. The scatter is also in- 
creased (0.53 dex, computed over the whole sample; 0.61 
and 0.55 dex in the small Mbh and small Mbuigo subsets 
respectively). We stress that the correction reported in 
equation [2] is uncertain, because of the wide range of 
B/T values reported for any given Ug. However, we re- 
mark that any correction for the B/T would make the 
case against a population of under-massive black holes in 
very massive galaxies even stronger. 

A similar argument can be used to evaluate how our 
results are affected by different assumptions on the star 
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TABLE 4 

Average values of the AfenZ-'^^host ratio in various subsets of our 
sample. (1) Subsample. (2) Number of objects. (3) Average value 
of log r. (4) Root Mean Square of log T. 



ii- The slope of the MBH~-Mhost relation in quasars is 
consistent with unity (in a log- log plane), consis- 
tently with what observed in quiescent galaxies. 



Su-bsctniplG 


N 
















(21 


(3) 




Without B/T correction: 








ALL 


87 


-2.842 


0.444 


AfBH > 10'' Mq 


33 


-2.627 


0.351 


Mbh < 10^ M0 


54 


-2.974 


0.445 


Mho.t > 4 X 10" M0 


49 


-2.889 


0.401 


Mho.t < 4 X 10" Mq 


38 


-2.782 


0.492 


With B/T correction: 








ALL 


87 


-2.697 


0.524 


Mbh > 109 M0 


33 


-2.623 


0.356 


A/bh < 109 Mq 


54 


-2.886 


0.603 


A/bulgc > 4 X 10" Mq 


46 


-2.886 


0.419 


AfbulKo < 4 X 10" Mq 


41 


-2.485 


0.554 



formation history. In our study, we adopted the mass- 
to- light ratio (M/ L) of a single stellar population origi- 
nated at Zburst = 5 and passively evolving down to z = 0. 
However, objects with significant disc contaminations are 
expected to have a younger stellar population than old, 
passive spheroids. This would imply smaller M/L, i.e., 
less massive host galaxies for a given observed host lu- 
minosity. This would make the case against a quasar 
population lying below the observed Men-Afhost relation 
even more robust. 

6. CONCLUSIONS 

We measured black hole masses and host galaxy lumi- 
nosities in a sample of 25 low-redshift {z < 0.5) quasars 
selected to have modest (< 10^ Mq) black hole masses. 
For each object we inferred stellar masses. This allowed 
us to significantly expand the sampled range of Mhost 
and A/bh for quasars. We found that: 

i- The MBH-Afhost relation holds over all the 2 dex 
both in terms of Mbh and Mhost- The relation has 
a scatter of 0.44 dex, i.e., the sampled parameter 
space is ~5 times larger. 
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iv- After applying a simplistic correction for the disc 
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V- No evidence of a population of quasars with rela- 
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found. 

Further studies at even lower Mbh masses (< 10^ M©) 
and at higher redshift could provide further constraints 
on the early black hole growth and the nature of the 
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galaxy relations along the Cosmic Time. This requires 
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that would become possible with the next generation of 
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Fig. 4 . — The AfBH^Afhost relation in the objects in our sample (filled symbols), as compared with the z < 0.5 sample from I Decarli et aLl 
H2010b"l (empty symbols) and the local relation (MBH/Afhost=0.0015; see Marconi & Hunt 2003). Squares, triangles and circles refer to 
A^BH estimates derived from H/3, Mgll and Civ respectively. Our new data substantially extend the sampled range of the MBH^Afhost 
relation. The relation holds down to quasar host mass es of ~ IO^^Mq (BH masses of ~ 10* Mq). There is no evidence of a population 
of objects lying below the relation, as claimed bv Batcheldoil II2010I) . The overall scatter of the MBll-M^^^^^ relation is 0.44 dex, and it 
extends over 2 dex both in terms of host galaxy and BH mass. Dotted, horizontal lines show how our objects would move if the correction 
for the bulge-to-total luminosity ratio is taken into account (see the text for details). The solid line is the best fit to our data. 
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